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2.1G FRESNEL ZONE CONSIDERATIONS FOR REFLECTION AND SCATTER 
FROM REFRACTIVE INDEX IRREGULARITIES 

R. j, Doviak and D. S. Zrnic* 

National Severe Storms Laboratory 
Norman, OK 73069 


INTRODUCTION 

Several different echoing mechanisms have been proposed to explain VHF/UHF 
scatter from clear air. GAGE and BALSLEY (1980) suggest three: (1) anisotropic 
scatter; (2) Fresnel reflection, and (3) Fresnel scatter, in order to account 
for the spatial (angle and range) and temporal dependence of the echoes. 

ROTTGER (1980) proposes the term "diffuse reflection" to describe the echoing 
mechanism when both scatter and reflection coexist. We present a unifying 
formulation incorporating a statistical approach that embraces all the above 
mechanisms and gives conditions under which reflection or scatter dominates. 
Furthermore, we distinguish between Fraunhofer and Fresnel scatter and present a 
criterion under which Fresnel scatter is important. 

Scatter from anisotropic irregularities of refractive index n has, for many 
years, been thought to be principally responsible for microwave echoes from the 
clear air. Existing formulations assume that the correlation length of n 
irregularities generated by turbulence, are small compared to the Fresnel 
length. But there is experimental evidence that the contrary may be true. This 
paper extends the existing formulations for the case where the Fresnel zone 
radius is comparable to or smaller than the correlation length. LIU and YEH 
(1980) recognized the limitations of the existing formulations which are based 
upon first-order expansion of the phase term in the integral for the scattered 
(or reflected) electric field intensity and WAKASUGI (1981) suggested expansion 
to second order. 

THE FRESNEL TERM IN THE INTEGRAL FOR ECHO POWER 

TATARSKI (1961, sect. 4.2) derived a formula for the field scattered from a 
volume V s with dimensions that implied the size of V s cannot be determined by 
the radar f s resolution volume Vg (DOVIAK and ZRNIC* , 1983). (The subscript 6 
is used to denote a resolution volume circumscribed by the surface giving a 
weight, to the scatterers, 6 dB less than the peak at the volume origin; DOVIAK 
et al • , 1979.) In a later publication, TATARSKII (1971, sect. 2.8) extended his 
earlier formulation so that V s could equal Vft. It can be shown that for back- 
scatter the condition assumed in this extension is 

P t « Ar s / 2ir - f/ (1) 

where r s is the range to an element of the scatter volume A, the radar wave- 
length, f the first Fresnel zone radius, and p t is the correlation length of 
refractive index irregularities for lags transverse to r s . Inequality (1) 
imposes the condition that constant phase surfaces of the incident wave are 
planes over the distance Pt> and the receiver is in the far field of this 
correlation length (i.e., r s > 2p t 2 /X). 

We now develop the scatter equations which allow correlation length to be 
larger than that specified by (1). Assuming the Born approximation (i.e., 
single scatter theory), the field intensity E-^, backscattered by refractive 
index irregularities An in the antenna far field is 
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where is the range to An at r (see Figure 1), f 0 (r) is the angular pattern 
of the incident electric field intensity assumed to be circularly symmetric 
about the beam axis, n 0 is the free space wave impedance, 1*^=2^/ A, P fc is the 
transmitted peak power, g the antenna gain, and r the distance from the origin 
of to the scattering element. V s is a spherical shel^ of thickness cx/2 
where t is the transmitted pulse width and therefore E^r^t) is the intensity 
of echoes sampled at a range-time delay (2r Q /c)+T after the transmitted pulse. 
Assume An to be a zero mean random variable. In a matched filter receiver 
having an internal resistance R, the increment of current magnitude |dlj 
produced by the scattering element is 
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where W(r) is the range weighting function (ZRHIC' and DOVIAK, 1978; DOVIAK and 
ZBNIC' , 1979). The integration now extends over all ? for which Wfg An has 
significant value. For a receiver filter matched to a rectangular transmitted 
pulse, the range weighting function is 
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where a is the unit vector from the origin of V r to the radar, 
o o 

The received power, time averaged over a cycle of the transmitted 
frequency, is: 

F = i II*R (5) 

r 2 

where * denotes the conjugate. 


From the integral of (3) 
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For the condition c T « r Q , r s does not change significantly where W(?) is 
appreciable so r g in the denominator of (6) can be replaced with r Q . Upon 
substituting (6) into (5) and taking the ensemble average, the expected received 
power becomes : 


<P r > * -4-4j{ R (r,r') W(r) W(r')f e 2 (r)f 0 2 (r' ) e" j2k 0 (r s _r s' } dVdV' 
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Let's assume th^t_^the irregularities have homogeneous statistical properties so 

that R(r,r* )-R(r-r f ) and that the two-way pattern function f 2 (r) is given by 

0 

£q 2 (?) “ exp {-e 2 M 0 Q 2 } 


(9) 
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where o Q 2 is the second central moment of the two-way jower pattern and 0 is the 
angular displacement, measured at the radar site, of r from the origin of V^, 

In terms of the 6 dB angular width 85 for the two way pattern f _ » 6^=3. 33a Q . 
For the assumption of narrow beams (9) can be approximated c 

f 2 =s exp {-t 2 /4r 2 a 2 ) (10) 

8 o v 

/ 5 2 

where t=/t + t ' is the projection of r onto the plane transverse to the beam 

axis at r 7 

o 

The Gaussian matched filter provides the best resolution of all the 
receivers having the same bandwidth (ZRNIC* and DOVIAK, 1978) . Because of this 
and because practical “matched filters" used in Doppler weather radars are 
Gaussian we assume that W(r) is well approximated by 

-K 2 i. 2 

W(f) « r ( 11 ) 

where is the second central moment of the weighting function W (r) and 

a ■« 0.35 cx/2 = 0.30 r, ( 12 ) 

r o 

for a Gaussian filter hatched" to a rectangular pulse of width t. The 6 dB 
range resolution is r^. 

Use the Taylor expansion for r g 

r .= r -r » r -a -r + {r -(a *r) } (13) 

s o 1 o o z r o 

o 

for terms up to second order in r 2 . We note that the second term of this 
expansion is the projection of r onto the a direction and the third term 
contains the projection on the transverse plane. Thus, in terms of the £,t 

coordinates centered in V , 

6 

r = r + t + t 2 / 2 r (14) 

SO o 

This quadratic expansion is valid (i.e., third order terms in r are negligible) 
provided that the scatter volume V g size is by limited by 

d t 2< 2Tq \/f 2 /n+d f 2 - |dj} (15) 


where 2d fc and 2d^ are the dimensions of V g transverse and parallel to r Q . The 
condition (15) assumes d^«r . As can be deduced, the farther the integration 
variable is displaced from the plane £=0 , the smaller must be the scatter 
volume size perpendicular to the beam axis. However, d^ in (15) need not be 
larger than the smaller of t^/2 or the longitudinal projection 
(d /cosii)+r Q 0 6tan$y2 of the scattering layer within (Figure 1). Substituting 
(14), the integral in (7) becomes 
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where a =a r Q /2 is proportional to the arc length of V . We now still find it 
convenient 0 to define new coordinates: b 
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Figure 1. Geometry for backscatter. The distances i and t are 
measured from the origin 0 of the resolution volume in directions 
parallel (longitudinal) and perpendicular (transverse) to the 
beam axis r Q . t^ is parallel to the x axis and perpendicular to 
t 2* 
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so that the tp t^ 1 component of (16) 


can be written as 


where <$ - 
from (16) 
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a l 5 l + ^2 <s 2 +a 3 5 3 * s tlie lag vector (note a^*^). The transformation 
to (18) is valid if, as is assumed here, the limits of integration 
entire volume where the integrand has significant value. Thus 
the integration over 


I 7 (t l ) " a t /iF 


R it) exp {-k Q ^ o t /r o ) 2 - 6^/40 
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Applying similar procedures to the t^ and l coordinate integrations we obtain 
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Resolution Volume Fresnel 

2 2 2 Weight Term 

where 6^= 6- + 6^ . The solution (20) is acceptable if the 2nd order expansion 
of r s in (147 is valid. Inequality (15) is the condition on V s for this 
expansion to be applicable. However, when the transverse dimension of V s is 
large such that (15) is not obeyed, we can still use (20) if R( o) is small when 
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~6 has a magnitude comparable to or larger than the right side of (15). In 
other words the correlation length p^_ perpendicular to the beam axis must be 

p t 2 < 2r Q |dj} (21) 

If condition (15) is satisfied, then there is no condition on p fc . By comparing 
(21) with (1), it becomes evident that the second order expansion relaxes the 
limits placed on the scatter volume size and correlation length p fc . Now these 
limits are increased by the factor (BttTq/ F or example, if r Q ~10 km and 
A=6 m, p t would have to be less than 1 .4 km in order for (20) to be applicable 
whereas would have to be less than 100 m for a 1st order theory. 

In the integral (20) the correlation is multiplied by two exponential 
weighting functions: (1) the resolution volume weight which depends solely upon 
the width a and range resolution cr r of and (2) the Fresnel terms which 
gives a weight in the t direction that depends upon t he ra tio f/ a t # Only when 
the radius of the Fresnel zone is large compared to /rra t p t can the Fresnel term 
in (20) be ignored. Therefore, both beam width and correlation length enter into 
the comparison with f . But because cr^ is a function of f , that is 


a 


t 


0.45 r 

< 

D/in2 


2 

0.9f 


( 22 ) 


where 0 is the antenna diameter, we can simplify the conditions so that the 
Fresnel term can only be ignored if satisfies 

p < (23) 

.t 0.97T 

On the other hand, because f is always smaller than a fc in the antenna's far 
field, the Fresnel term in (20) will have more weight than the beam width part 
of the resolution volume term. Thus situations that allow us to neglect the 
Fresnel term will also permit us to ignore beam width influence. If (23) is 
satisfied, we can use (20) (without the beam width and Fresnel terms) to obtain 
the scattered field, even though V s is larger than V^; then we need to sum 
incoherent echo power from elemental volumes large compared to p ^ but small 
compared to Vg (DOVIAK and ZRNIC* , 1983).' We call this case incoherent 
Fraunhofer scatter. But H0DARA (1966) shows that within the lower troposphere, 
the correlation length has the following height dependence 


p * 0.4h/(l+0.01h) (m) (24) 

Where h is in meters. Furthermore, VHF backs cat ter data analyzed later in 
this paper suggest that p t ~20 m for irregularities in the lower stratosphere. 
Thus, unless the antenna diameter is of the order of 100 m or more, the Fresnel 
term will be important in determining the field scattered by refractive 
irregularities. If the scattering volume contains many subvolumes for which 
(20) applies, but (23) is not sa tisfied, we have a situation of incoherent 
Fresnel scatter. When d t </2r Q f (from Equation 15), then we have coherent 
Fresnel scatter. If d <f , then signal is coherent irrespective of the trans- 
verse reshuffling of refractive index irregularities. 


THE SPECTRAL SAMPLING FUNCTION 


Because it is common to describe the statistical properties of refractive 
index irregularities by the spectral density function the effects of the 
resolution volume and Fresnel terms on echo power can be examined conveniently 
by introducing a spectral sampling function. Equation (20) can be expressed in 
terms of the Fourier transform of R(6) multiplied by the lag weighting function 
H(l) where 
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Thus 
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I 7 = $ nw <0 '°» 2k o ) (26) 

where $ nw (K) is the three dimensional transform of R(6) multiplied by the lag 
weighting function. Now $ nw is the spectrum $ n of refractive index 
irregularities convolved with the spectrum of H( o) : 

$ = 8 0 0 2 it 9/ ’ 2 ft $ * $ (27) 

nw r t n w 


where * denotes convolution and 
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is the normalized spectral sampling function. Substituting (25) into (28) and 
evaluating: 
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where K t = /kT^+K^ . The second order phase term has contributed the factor 
4u f 4 in 1 the above equation. Thus the first order expansion is valid only 

if this factor is small relative to unity. However for in the antenna far 
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For remote sensing with radar it is^commonn to have in the antenna far 
field, thus the Fresnel term in $ W (K) cannot be ignored. As discussed earlier, 
this conclusion is a result of the fact that the Fresnel radius is always less 
than the beam width so that the Fresnel term always dominates the beam width 

weighting function. Thus $ (it) can be well approximated by 
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in which we have substituted (22) for a t . Equation (31) shows that the larger 
is the antenna diameter, the narrower is the spectral sampling function. It is 
surprising that the sampling function shape and size is independent of r c and, 
for a given antenna diameter, the spectrum $ n (ft) of irregularities is weighted 
equally for all resolution volumes in space. This result differs from that 
derived by TATARSKII (1971) who only considered first-order phase expansion in 
which case is a function of r Q . By combining (7), (26) and (27) the back- 
scattered power is given by 


<P > 
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In the atmosphere it is usual for the horizontal correlation length p n to 
be larger than the vertical one p z so $ (ft) will be more sharply peaked along 
the K x ,K y directions and less so along the K z axis. If the irregularities have 
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shapes that are roughly described as oblate spheroids, then the correlation 
R(^) would also have a similar form but $ (it) would be prolate spheriodal in 
shape (Figure 2a). Equation (31) reveal s n that whenever range resolution 
rg=3.33a r is larger than 0.34D, as is usual, the sampling function (Figure 2b) 
along will be narrower than along K t * If the beam axis is rotated by ^ 
degrees from the vertical, $ w (£) will also be rotated by $ from the K z axis. 

If echo power decreases significantly as iji is increased, then we have 
specular type reflection. The sharpness of the angular dependence is a function 
both of p and D. Referring to Figures 2a, b and Equation (31), we see that a 
necessary n condition to observe a specular type echoing mechanism is for 
0.54\/D«l. That is, narrow beams are required which is consistent with simple 
physical arguments. Assuming space is filled with An, specular type echoes 
will then be observed only if p^»P • However, we must be cautious in applying 
these criteria because we have used a specific model (i.e., Gaussian) to 


(a) (b) 




Figure 2. (a) Contour surface of constant spectra intensity $ n (K) 
for irregularities having symmetric correlation lengths along 
x and y that are longer than the correlation length along z. 

(b) Contour surface of the spectral sampling function $ W (K) for 
beam axes at^elevation angle G e = tt/2 - ty. (c) Contours of ' 
constant $ n (K) for which the small-scale irregularities produce 
isotropic scatter. 
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describe the statistical properties of n and because Ph,P z only characterize 
the most intense irregularities of refractive index* Thus, if the contours of 
constant $ n (*b have the dependence sketched in Figure 2e, scatter could be 
independent of ^ (i.e., isotropic scatter) if 2k 0 »p z ~^ and 2k 0 »v5r^ -1 . The 
$ (it) depicted in Figure 2c can be represented by a sum of isotropic <5 . and 
anisotropic $ a p^rts whei^e is the wave number beyond which (independent 

of direction of K). $..(K) could have the -11/3 power law dependence on K 

deduced from turbulence theories. 

BACKS CATTERING FROM ANISOTROPIC IRREGULARITIES 

As an example, let us consider the angular dependence of echo power when 
the scattering medium can be decomposed into isotropic and anisotropic 
components (i.e., R(6)=R i +R a ). We further assume that R is isotropic in the 
horizontal plane. To obtain an order of magnitude estimate, we take R of the 
form: ^ 2 ^ 2 d 

R = <An 2 > exp {■ — — 77 - • Z ) 


6 , = 6 2 (33) 

h x y v ' 

The resolution volume coordinates are related to the natural coordinates x,y,z 


6 =6, ; 6 “6 cos^-6 0 sinf; 6 sinip+6 0 cosijj (34) 

x 1 y 2 3 z 2 3 

After introducing (34) into (33) and the result into (20), integration is 
performed giving the formula for echo power from anisotropic irregularities as 
being proportional to 
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we can simplify (35) 
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(39c) 


When height of V, is constant and laminae are infinitesimally thin, the ratio 
for echo powers at angle ip and zenith ('HO) is: 


<P a («(')> 

<P a (o) > 
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The term cos $ accounts for the decrease in power due to the range and 0 t 
increase with tilt away from the vertical because remains at constant 
height. To <P we add the power due to isotropic irregularities to 

obtain: a 
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<p(n> 
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where A-P i (0)/P a (0) ♦ Equation (41) was fitted to data (Figure 3) from ROTTGER 
et al, (1981). Pertinent parameters for the Rottger et al. data are: X=6.4 m; 

D=*260 m; heights h=16.9-*18.1 km near the tropopause ; beam width 01=1.7°; and 
range resolution = 300 m. We find P h =20 m for the horizontal correlation 
lengths, and A=0.04 fits well these data. Comparing terms in (39) it is seen 
that the Fresnel term (i.e., the 2nd term in Q(^)) does not contribute 
significantly. Although we have not distinguished any one of the mechanisms 
discussed in the introduction as being responsible for the echo power, we see 
that scattering from anisotropic irregularities can account for the observed 
angular dependence which is sufficiently peaked that one might believe a 
reflection mechanism is acting. 


For sake of simplicity, it is preferable to label the echoing mechanism as 
scatter whenever there are several or more scattering irregularities for which 
only a statistical description of their properties (e.g., size, intensity, etc.) 
is practical. Thus, we do not need to invoke a reflective process to explain 
observations in this case; the scatter formulation presented here can explain 
all the features of the received field if indeed the medium is comprised of many 
irregularities of refractive index for which only statistical properties are 
known. 


Equation (2) is the starting point for our formation for scatter from 
refractive irregularities. Although we refer to (2) as the scatter integral, it 
can be used as well in situations (i.e., P t »f) which might be interpreted as 
reflective. In order to determine echo power when irregularities have 
horizontal dimensions large compared to the Fresnel radius, GAGE et al. (1981) 
have used the general formula 
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TILT FROM VERTICAL ir/2-9 (deg) 


Figure 3. Angular dependence of observed mean backscatter power 
(open circles) from anisotropic irregularities as the radar 
beam axis is tilted away from the vertical (ROTTGER, 1981) . 
Fitted to the data is a model that consists of anisotropic tur- 
bulence with a two-dimensional (horizontal isotropy) correlation 
function in an isotropic background. 


1/2 

1* I 2 “ 4 i n di eXp { “j 2k 0 z } dz 1 2 (42) 

-i/2 

for the power reflection coefficient where l is the thickness of the partially 
reflecting layer* In this form variations of n along the horizontal are ignored 
and, if the scattering layer is in the antenna far field, the echo power P r is 
easily found by considering an image source which gives 

P 'PA 2 |R| 2 /4A 2 r 2 (43) 

r t e 1 1 o 

where A e is the effective area of the antenna (A e -gA^/4). For exactly the same 
assumptions on n, the solution of (2) should produce an identical echo power. 

In Appendix A we prove this contention by simply using the second-order phase 
terms; this shows the wide applicability of the solution presented earlier. 

Figure 4 illustrates the type of scatter that would be effective versus the 
location of the sampling wave number 2k 0 for the case P^»P Z « The location of 
boundaries are functions of the parameters p h , D, r Q and the relative strengths 
of $ a and and thus there could be a different order than presented on 
Figure 4. For example, if K i <2r 0 /p h , then Fraunhofer scatter could be either 
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Figure 4. Types of echo mechanisms versus the location of the 

sampling wave number 2kQ for a particular ordering of boundaries. 
P\ l »Pz> an <i v 6 is uniformly filled with irregularities. 


anisotropic or isotropic depending upon the value of 2 k * Because we had 
assumed infinitesimally thin laminae (K /* 00 in that case), scatter will be 
anisotropic no matter how large is 2k Q . However, it is more likely that 
K i < 2 'frr 0 /p h so that if 2 k 0 was larger than 2 m " 1 we could pass into a region of 
isotropic scatter. Data collection at various 2k o * s (i.e., multiple wavelength 
radar) could establish the correlation length p z _^and a value for K,. At UHF 
wavelengths 2 k 0 is so large^that the peak of $ w (K) is expected to fall most of 
the time in the tail of $ n (K) at wave numbers where turbulence is mostly 
isotropic and <f> n is expected to have the same ll/3rds dependence on K as does 
the velocity fluctuations. However, at the longer wavelengths in the VHF band, 
2 k 0 is much smaller so it can place the $ w (it) peak in a region where may 
sometimes be larger than ^ or smaller than it. 

ECHO POWER DEPENDENCE ON RANGE AND RANGE RESOLUTION 

GAGE et al. (1981) propose a model for which echo intensity varies as the 
inverse square power of range but has a range resolution dependence that can 
vary from zero ‘to a square law. BALSLEY and GAGE (1981) introduce the concept 
of a scatter volume defined, transverse to the antenna beam, by a correlation 
radius to derive an echo intensity that depends on the fourth power of range. 

It is improper to form such a condition because the scatter volume V s is 
defined by either the spatial distribution of intensity of An fluctuations or 
by the resolution volume V 5 , whichever is smaller. We shall use the solutions 
derived here to determine the conditions under which various dependences can 
occur. Recently HOCKING and ROTTGER (1983) have critically reviewed the inter- 
pretations of Balsley and Gage. 

Assume vertical incidence and use (31) and (32) to obtain 

Co _y ..y 

<P > = — | 0 (K)* $ (K) (44) 

r 2 n w 

r 

o 

where C is a constant independent of and r Q , and $^(1?) can. be expressed as: 

$ ($) - o $ <K ) exp {-2a 2 K 2 } 
w r w t r z 


( 45 ) 
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2 2 2 zt 

where now K « K + K , Consider two cases: (1) $ (k) broad and (2) narrow 

t x y n 

compared to $ (it) along the wave number it coordinate (Figure 5). 

w z 

(a) $ Broad 

n 

Integration along gives a <P^> approximated by 


$ (K ) $ (K ,K ,2k )dK dK 
w .t n x* y* o x y 


(46) 


which illustrates that the expected echo power is proportional to range resolu- 
tion (assuming a uniformly filled V^) and inversely proportional to the square 
of range r Q . This ia the usual dependence expected when scatter is from 
irregularities produced by turbulence. The r “'dependence occurs irrespective 
of whether is large or small compared to f ? 

(b) $ (it) Narrow 

n 


In this case (44) can be reduced to 

c„ a -vw 2 

<P > * e 

r 2 

r 
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* (K ) $ (K) dV^ 
wen k 


Again <P r > depends upon the inverse square of r G ,a result which is independent 
of p^* If $ w (K t ) is also broad compared to $ n C&) along K t then (47) reduces to 

C0 2 -20 2 (K -2k ) 2 r 

<P r > = — f~ - V 0)e r SZ ° j $ n (K)dV K (48) 

r o 

in which $ n (lt) is assumed to have a peak at K *0* This case occurs when 
refractive index irregularities have pft>f and strong Fourier components 
clustered about 2k Q . Only if 2k 0 =k sz will <Pj> be proportional to other- 

wise, we could have other range resolution dependencies. The integral in (48) 
is the variance <&i 2 >. 



Figure 5. Cases in which $ n is broad (a) and narrow (b) compared to 
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If the irregularities are contained in a thin layer having a vertical 
dimension small compared to a r , then echo power would be independent of tf r , 
when the resolution volume is centered on the layer. For resolution volumes 
displaced from the layer height, echo power would have a strong dependence on 
a r even exceeding the square law one! If the scattering irregularities are 
confined to horizontal dimensions small compared to beam width, then a fourth 
power range dependence would be obtained. However, the correlation length does 
not determine the scatter volume dimension as stated by BALSLEY and GAGE (1981). 
In the cases discussed in (a) and (b) the echo power depends on the inverse 
square of range because we have assumed uniformly filled V6 . 

CONCLUSIONS 

When scattering layers are in the far field of an antenna, the Fresnel term 
is a more important weighting function than the antenna pattern because the 
width of the antenna pattern is always larger than the Fresnel radius f. Only 
in the case where the correlation length p of refractive index irregularities 
to perpendicular to the beam is much smaller than f will the first-order 
truncation of the Taylor series expansion for phase be valid. Then Fraunhofer 
scatter is considered to be effective. However, when retention of second-order 
phase terms is necessary a Fresnel term (see Equation 20) is introduced. The 
criterion for keeping the second-order phase term depends both upon beam width 
and the Fresnel radius. Thus, the condition under which incoherent Fraunhofer 
scatter is effective becomes solely a function of antenna diameter D (i.e., 
P t <0.29D). When the Fresnel term needs to be included in the solution we have 
the situation of Fresnel scatter or reflection. It is suggested that unless the 
antenna diameter is of the order of 100 m or more, the Fresnel term is important 
in determining the field scattered by refractive irregularities. 

The formulas derived here establish the conditions under which a scatter or 
reflection mechanism can be distinguished. However it is important to have the 
proper statistical description of the irregularities in order to obtain the 
spatial and temporal dependence of echo intensity. A multiple wavelength radar, 
in which its beam position can be scanned, could supply invaluable data to 
characterize the spectrum of refractive index irregularities and help to explain 
the properties of the echoes. Only when irregularities have a spatial spectrum 
f orm that concentrates variance <An^> at wave numbers near 2k 0 does echo power 
depend upon the square of range resolution. Echo power depends upon the inverse 
square of range r Q independent of whether p t is less than or greater than f. 
However the resolution volume must be uniformly filled with An. 
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APPENDIX A 


In this appendix we shall demonstrate that (42) and (2) and identical 
formulations for the situation considered in this paper (i.e., a scattering 
layer in the far field of an antenna, and range resolution sufficiently narrow 
so that the 1/r^ term can be brought out of the integral (2)). 


For pulsed transmissions the height interval that contributes to the echo 
sample is determined by the pulse shape if refractive index Irregularities are 
distributed throughout the vertical. Therefore, in this case (dn/dz)/n in (42) 
must be multiplied by the range weighting function (e.g. , (11) ) and then, for 
pulse widths small compared to r Q , the limits on z in (42) can be increased to 
infinity without significant error. Thus the reflection coefficient takes the 
form: 



-z /4cr 

r 
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— (£nn) exp {-j2k z} dz 
dz o 


—oo 


(Al) 


Using integration by parts and noting that n - 1 + An where An«l, (Al) can be 
reduced to 
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2 

Now for range resolution many wavelengths long (i.e., k 0 o r »l) the term z/4a r 
can be ignored in the integral without adding appreciable error to p. Then 


R - jk 
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oo 

» 

W(z) An(z) exp {~j2k z } dz 
0 

—00 


(A3) 


We note here that W(z), as defined in this paper, also contains the weight 
associated with the frequency transfer function of the receiver's filter. 
Although this function does not rigorously belong in the integral for p, we 
are primarily interested in the received echo power which is dependent upon the 
filter function. For a linear system we could, if we ignore receiver noise, 
just as well consider the filter at the transmitted output thus modifying the 
pulse shape to give the equivalent weight W(z) considered herein. With similar 
consideration we can also express (2) in the form 


^r)^ exp {-J2k r } dV (A4) 

o s 

We now consider An to depend upon z as in the reflection formula and using (10) 
for f 2 (r), the second-order expansion (14) for r , and integrating over the 
horizontal we obtain ° 
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W An exp {-j2k o z}dz (A5) 


Now for r Q in the antenna far f ield the term jk^r Q has a magnitude larger than 
1/2 a t 2. The echo power P^ in terms of E ^ is: 

P r “ A el E l| 2/2n d (A6) 


substitution of (A5) and (A3) into (A6) reveals 


P 

r 



which is identical to (42) derived from the reflection formula. 



